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As a first step toward circumventing the difficulty to obtain an attractive isospin-independent 
AW spin-orbit force from Skyrme-type models involving only pions, we investigate an improved 
Skyrme Lagrangian that incorporates the scalar-isoscalar meson e which can be viewed as the cause 
behind the enhancement of the 7T7t S-wave. We find that at large distances, the main contribution 
to the spin-orbit potential comes from the scalar Lagrangian and it is found to be attractive. We 
^\ ' briefly discuss how to pursue this work to finally obtain a medium-range attractive interaction. 
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I. INTRODUCTION 



Ten years before the advent of QCD, Skyrme proposed a model Q| describing hadronic physics which involves 
only pion fields and where baryons emerge as topological solitons. This model is recognized as the simplest chiral 
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realization of QCD at low energies and large N c The corresponding Lagrangian contains, in addition to the 
i- £h ' well known nonlinear a model, an antisymmetric term of fourth order in powers of the derivatives of the pion field 
(the so-called Skyrme term). The latter has been added by Skyrme in order to avoid soliton collapse. Despite its 
relative success in describing some properties of baryons (cf. Ref. [0 for a review), the Skyrme model presents several 
shortcomings such as its prediction of a very small value for the axial coupling constant (J] and an almost zero nucleon 
mass when the Casimir effects are taken into account ||. But probably the most important crisis of the Skyrme 
model is that it does not allow for the formation of nuclei. Indeed, it predicts a repulsive NN central potential ||. 
Because of these drawbacks, several extensions of the standard Skyrme model have been proposed in order to improve 
its prediction power. The most famous one consists in adding to the Skyrme Lagrangian a term of order six in powers 
of the derivatives of the pion field, proportional to the square of the baryon current as first proposed in Q). The 
extended Skyrme model thus obtained contains a second, fourth and sixth-order term, but still is an effective theory 
of pions. It has been shown recently [j5],^| that this extended model is much more realistic when describing low-energy 
hadronic physics than the standard Skyrme model. 

In contrast with the velocity-independent part of the NN potential, the kinetic part and especially the NN spin- 
orbit force has not been studied extensively within the framework of Skyrme-type models. However, all the authors 
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who worked on this subject, e.g. |s|JlC(l, arrived at the disappointing result that the standard Skyrme model predicts 
a repulsive isospin- independent spin-orbit force while, according to the phenomenology |p"l|| , it should be attractive. 
A glimpse of hope has emerged from Ref. fl2]| where it was claimed that when considering the sixth-order term, an 
attractive isoscalar spin-orbit interaction could be obtained. However, we have shown recently [ fL3| that this is not 
the case: the sixth-order term contributes with a positive sign as is the case for the Skyrme term. Indeed, the authors 
of Ref. jl2| considered in their calculations only one piece of the force due to the sixth-order term and omitted the 
piece that stems from the baryon exchange current, which turns out to be not only repulsive, but the dominant one 
too jn3. Thus, despite its success in improving the predictions of 1-baryon data, even an extended Skyrme model 
including higher-order terms in powers of the derivatives of the pion field fails to predict an attractive NN isoscalar 
spin-orbit interaction as it fails to predict the correct NN central potential Therefore, one has to consider more 
realistic low-energy models to describe the NN interaction. 

The Skyrme term as well as the sixth-order term can both be regarded as an infinite mass limit, the so-called local 
approximation, of a model with p-meson [ [l5[ and w-meson H exchanges respectively. These findings have provided 
the ground for improving Skyrme-type models by constructing effective Lagrangians which include, in addition to 
the pions, the other low-lying mesons. The authors of Ref. have investigated the NN force by considering such 
an effective model incorporating the mesons 7r,p, lu,Ai and the scalar meson e. Unlike previous works, they have 
introduced the vector meson fields by gauging the linear er-model instead of the nonlinear one. They have found 
that the different channels of the NN interaction arising from the potential piece of the Lagrangian are very well 
described, in particular the central potential for which an attraction has been found in a quantitative agreement with 
the phenomenology These results lead us to think that, similarly to the central potential, it may be that the 

right way to circumvent the difficulty of finding an attractive isoscalar spin-orbit force with Skyrme-type models is to 
consider effective theories which incorporate low-mass mesons with finite masses. In order to check these claims, one 
has to investigate the kinetic energy of a two-nucleon system within an effective Lagrangian such as the one in [fl6[ . 
This is of course not an easy task and goes beyond the scope of this short paper. As a first step, we rather focus here 
only on the scalar meson contribution to the isoscalar spin-orbit force. Indeed, it has been pointed out in Ref. fjlf l 
that the presence of scalar degrees of freedom plays the crucial role of recovering the missed central attraction, so one 
may think that this will also be the case for the isoscalar spin-orbit force. 

This idea to account for a scalar field in the spin-orbit channel has been investigated in Ref. p7|| . There, a dilaton 
field is coupled to Skyrmions in order to mimic the scale breaking of QCD. We will not discuss here if a dilaton 
field is suitable or not to provide a good description of low-energy hadron physics; we refer the interested reader to 
Ref. [|l8| for a detailed discussion on this subject. The fact is that a combination of the sixth-order term and dilaton 
coupling yields an attractive isoscalar spin-orbit force, as it has been claimed in p7| . Note though that this result 
remains questionable since these authors, like those of Ref. [fllf , ignored the baryon exchange current contribution in 
their calculations |jl3|| . In any case, here we will follow a more transparent way in order to incorporate scalar degrees 
of freedom |jl6| . It consists in considering the linear a model in which the presence of the scalar field is intrinsic. 
In Sec. II, the model with the scalar-isoscalar e-meson is presented. Sec. Ill contains the derivation of the isoscalar 
spin-orbit force. Our results and discussions are presented in Sec. IV. 
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II. THE MODEL WITH SCALAR MESONS 



The infinite-mass limit of the effective Lagrangian with ir, p, u> and e mesons jlfj reads: 



C — £nlo- + C± + Cq + £ x sb 

f 2 1 9 f 2 m 2 C 1 ) 

= J -fTr (d li Ud»U+) + ^Tr [(d li U)U+, (d„U)U + ] 2 - e 2 B^ + ^_JLTr(U - 1) 



where = e^^Tr ( (d u U)U + (d a U)U+ {d () U)U + ) /24tt 2 is the baryon current §, and U an SU(2) matrix which 
characterizes the pion field. The first term in Eq. (|l|) corresponds to the nonlinear a model, f n being the pion decay 
constant. The second term, parameterized by the coupling constant e, is the so-called Skyrme term. The third term is 
of order six in the derivatives of the pion field and corresponds to w-meson exchange in the case of an infinite w-meson 
mass [Q. The coupling constant e\ is a parameter related to the u — > ttj width. The last term in Eq. (|l|) which is 
proportional to the square of the pion mass (139 MeV) implements a small explicit breaking of chiral symmetry. 
As compared to the standard Skyrme model, the extended model (Q), when used with realistic parameters, provides 
a more accurate description of the 1-baryon properties (e.g., the nucleon quantum mass, the A-N mass splitting, the 
breathing mode energy the axial-vector coupling constant gA,---) However, as mentioned in the introduction, it 

does not describe properly neither the NN isospin-independent spin-orbit force [fl3| nor the central potential Jl4| . 

The scalar-isoscalar meson e is viewed as the responsible for the S'-wave attraction in the itit interaction and, 
therefore, should be considered as an essential ingredient in the low-energy hadronic phenomenology. It is incorporated 
in the Lagrangian (pi) by replacing the nonlinear a model £nlct with the linear one. Namely, 



£nl. - £l„ - ! W£ + L Tr (d ll Ud*U + ) - X(e - flf ( 2 I 



where £ is the isoscalar content of the quaternion field a + iT.n, a being the scalar chiral partner of the pion field. 
The scalar field e is defined as e(r) = /„■ — £(r). The coupling constant A is related to the e- meson mass through 
m 2 = SXf 2 . It is straightforward to check that by taking m e — > oo in Eq. (|2|), one recovers the nonlinear a model 
£nlct. Indeed, in that limit, A becomes infinite and so the £-field has to be frozen to its asymptotic value / w in order 
to keep the potential energy in Eq. ((H) finite. 

Therefore, for a realistic finite e-meson mass, the Lagrangian density (Q) is then replaced by the following one: 

C = Cha + £4 + £e + -CxSB (3) 

The one-soliton system is commonly solved by assuming the hedgehog ansatz for the pion field: 

U{r) = U H {r) = exp ( iT.rF{r) ) (4) 

where the r a 's are the Pauli matrices and the notation f means v/r. The static Euler-Lagrange equations for the 
chiral field F and the £-field corresponding to the Lagrangian (||) read 

1 el 



2?r 4 r 2 



K)" = (V 2 + 4r sin 2 (^) + 4A(£ 2 - / 2 ) ) r£ 



F ') =^ 2 + ^sin 2 (F) jsin(2F) + 
sin 2 ( J F) J sin(2F) F 2 + / 2 m 2 r 2 sin(F) , (5) 
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where primes indicate radial-coordinate differentiation. In order to ensure a winding number one the chiral angle F 
obeys the usual boundary conditions while the £-field fulfills the conditions £'(0) = and £(oo) = f w . Let us observe 
that the above Euler-Lagrange equations are solved with a set of parameters whose values have been fixed by fitting 
to the mesonic sector |5|,|8|,|l6| . Namely, /„• = 93 MeV, e = 7.2, ee — 1.66 fm, m e = 650 MeV. These parameters yield, 
e.g., a value of the nucleon mass of ~ 1 GeV after subtracting the Casimir energy ||, and a A-N mass splitting 
of ~ 267 MeV. We have plotted in Fig.l the chiral function F and the £-field, solutions to the equations of motion (^) 
obtained with the set of parameters given above. We will use below these two functions as an input in the numerical 
computation of the spin-orbit force. For reference, we have also shown in Fig.l the chiral solution F in the case of 
the model (@, i.e. £ = /„,§. 

III. THE SPIN-ORBIT FORCE 

So far, all the calculations on the NN spin-orbit force within Skyrme-type models have been carried out in the 
framework of the product ansatz as suggested by Skyrme jl9| (cf. Ref. ||,|l0| and references therein). It is worthwhile 
noticing that the latter is only a simple approximation to the two-baryon configuration and not a self-consistent 
solution. We will not discuss here the degree of validity of this approximation. It is commonly chosen because of 
its relative simplicity as compared to other two-baryon field configurations which can be found in the literature p0[ . 
Furthermore, it becomes exact for large NN separation^. Thus, following common practice, we use it here to describe 
a system of two interacting solitons. Furthermore, in order to obtain the appropriate spin and isospin structure, we 
also introduce rotational dynamics [Q. Hence, the field configuration of the two- nucleon system separated by a vector 
r reads: 

U(A 1 ,A 2 ,x,r) = U 1 U 2 = A 1 U H {ri)A+A 2 U„(r 2 )A+ , 

(6) 

ri=x — r/2 , r 2 =x + r/2 , 

where A\ and A 2 are SU(2) matrices, and Uh the hedgehog single soliton (||). For the scalar field, we use the 
configuration suggested in Ref. jfj) 

fr(x,r) = i£( ri )£(r a ) (7) 
J n 

which obviously is the form the most compatible with the product-ansatz approximation (|^) . To carry out a simulta- 
neous quantization of the relative motion of the two nucleons and the rotational motion, we need to treat r, A\ and 
A 2 as collective coordinates. Therefore, we make all these parameters (r, A\,A 2 ) time-dependent. 

The spin-orbit potential will emerge from a coupling between the relative motion and the spins of the two nucleons 
so that we have to calculate the kinetic energy corresponding to (j^). As reported in Ref. 10 13], generally, one has to 



treat with care the conversion from velocities to canonical momenta before identifying and extracting the spin-orbit 
potential. Indeed, one has to start from a Lagrangian formalism, take its "classical" kinetic energy and extract from it 
the mass matrix and then invert it properly in order to move to a Hamiltonian formalism. However, for a large relative 



1 The region of validity of the product ansatz corresponds to a relative distance r much larger than 1 fm. 
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distance r, the region of validity of the product ansatz, this procedure is equivalent to that of Refs. |g,|l^,|3| which 
we will use here. In this latter one starts, in the case (||), directly from 

Kur{A u A 2 ,r) = (-1) J d 3 x (^e P Tr(d Ud a U + ) + i(cW) 2 ) (8) 

where the minus sign in front of the integral is put in explicitly so that account is taken of the change in sign of the 
off-diagonal terms of the 2x2 mass matrix under inversion 0. Afterwards, one makes the usual identifications : 

nW i . s («) 

' "n = -^(-rA+A n )^— , n=l,2 , (9) 

where and s^™) are respectively the radial momentum and the spin of the n-th nucleon while M and A are 
respectively the mass and the moment of inertia of the single soliton. Inserting the product ansatze ([K?]) in the 
expression (||) gives 

= J d 3 x ^ 2 (r!)e 2 (r 2 ) (Ro a {U x ) + L 0a {U 2 ) f +■■■ (10) 

where we have omitted to write the second term in Eq. (||) since it does not contain angular velocities and thus will 
not contribute to the spin-orbit force. In Eq. (|l0|), i? 0a and L 0a are the time components of the right and left currents 
respectively: 



Roai^) = - l -Tr(r a U+ cWi) - D ab {Ai) \Tg ^ + 
L 0a (U 2 ) = - l -Tr(T a d U 2 U+) = D ab (A 2 ) ( 




(11) 



(n) in) 

In the above equations, the sum from 1 to 3 on repeated indices is understood. The tensors Tj 6 and R ab depend 
only on the position r n of the n-th nucleon, and their explicit expressions can be found in Eq. (11) of Ref. ]l3[ ] . The 
D a bS are the matrix elements of the 3x3 rotation matrix in the adjoint representation: D a b(A) — Tr (r a ATf } A + ) /2. 



To obtain the right-hand sides of expressions (11), we have used the quantization scheme (|9|). The next step now is 
to expand the square in the expression of ( |l0| ) after inserting expressions ( |TT| ) and making use of the D-matrix 
properties and the definitions of the tensors and R^' |l3|. Thus, by considering the isoscalar part^J and keeping 
only terms proportional to L.S, where S = + is the total spin and L = r x p the angular momentum (p being 
the relative momentum, i.e., p = p' 2 - 1 = — p^), we obtain: 

MA l r l (12) 

V L „(r) = -75- / d 3 x £ 2 ( ri )£ 2 (r 2 )- sin 2 ^)) f.f x . 

The total isoscalar spin-orbit potential of the model (||) is then: 

Vso,z=o = ^ (V L «(r) + V 4 (r) + V e (r)) L.S (13) 

where the linear a-model contribution Vl<t is given in Eq. (]l|). The rather lengthy expressions of the Skyrme- 
term contribution V4 and the sixth-order term contribution Vg can be found explicitly in Refs. [j^,|l0|,[l7| and |L3| |, 



Due to the projection theorem, the terms proportional to D(A^ A2) contribute to the isospin dependent force. 
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respectively, and do not need to be rewritten here. Note that the analytical expressions of V4 and Vg are not affected 
by the presence of the scalar field and still depend explicitly only on the chiral field F. This is because the £-field 
does not couple to £4 nor £g, as it can be seen from Eq. (^). But obviously their numerical values will change due to 
the dependence of F on the scalar degrees of freedom. 

IV. RESULTS AND DISCUSSION 

Before displaying our numerical results concerning Vl<t , let us first discuss the analytical properties of the terms 
Vi(r) and Vq{t) at large distances [cf. Eq. (fl3|)]. As already shown in Ref. [ fL3| , these potentials decrease as e~ 2m ^ r a 
large r. Indeed, one can see from Eq. (||) that at large distances the chiral angle F(r) falls off as e~ ni7,r /r (note that 
at large r, £ — > /„■). Then, following the approximation of Ref. which consists of treating at large distances the 
field of one soliton as constant in the presence of the other, one obtains from [p 10,1711 that V4 decreases as e~ 2m * r 



over a power of r. For the potential Vq we have shown in Ref. |13| ] that it is a sum of two terms. The first term (the 
one which has been considered in p^ , p7| ) decreases as e _3m,rr . This behaviour is expected since the sixth-order term 
Cq in Eq. ([!]) can be regarded as the local approximation of an effective model of ui meson coupled to three pions . 
However the asymptotic expression of the second part of V§ 7 which arises from the baryon exchange current p3[ , 
contains both e - 2m ^ r and e -3 ™ 1 "'" terms so that it behaves at large r as e ~ 2m " r . This means that the Lagrangian Cq 
generates, in addition to the expected three-pion exchange piece, a two-pion exchange one (coming from the exchange 
current) similar to that of the the fourth-order Skyrme term £4. Thus we see that for large distance both potentials 
V4 and Vq are of Yukawa type with range 2771^. In fact, the main motivation of including a scalar- isoscalar meson 
field by hand in Skyrme-type models [ ]l6| , |l7t is to mimic the two-pion exchange since the latter is not well addressed 
in these models with simple zero-mode quantization. 

A straightforward manner to calculate numerically the integrals giving the different contributions to the spin-orbit 
force consists in making the change in variable r' = x — r/2 and taking the TV TV separation vector r along the z' axis. 
E.g., with these changes, the expression of the linear a model potential Vlo- in Eq. ( |l2"| ) becomes: 

/■°o /■! 
- / At 1 r' 9.™ 2 (V(r'\\f: 2 (r'\ I 

ft 

where u = r.r' — z'/r'. In Fig. 2, we plot Vl<t, V4, Vq together with the total sum Vl ct + V4 + Vq as functions of the 
relative distance r between the two nucleons. As it was the case for the extended Skyrme model (|l|), V4 and Ve are still 
repulsive in the case of the model @ . On the other hand, while the nonlinear a model £nl<t has a zero contribution 
to the isoscalar spin-orbit forced], the linear a model contributes and with a negative sign as it can be seen from the 
behaviour of in Fig. 2. This result is in agreement with the phenomenology [jll]]. The total potential (full line 
in Fig. 2) is found to be repulsive for intermediate distances (r < 2fm) and attractive for large distances (r > 2fm). 
For short distances (not shown on Fig. 2), neither the product ansatz nor the effective model should be trusted since 
this region corresponds to processes involving perturbative QCD. The total potential is attractive at large distances 



V ha (r) = I dr' r'sm 2 0F(r')KV) / u (, 2 ( Vr' 2 + r 2 + 2rr'u) (14) 



because in this region the scalar meson contribution (12) is the dominant one as compared to V4 and Vq. This can 



3 Indeed, £nl<t contributes only to the isospin-dependent force 
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be checked straightforwardly by comparing the asymptotic behaviour of each expression in Eq (|l3|). We have also 
displayed in Fig. 2 the potential related to the linear a model and the sixth-order term only, Vl<t + Vq. We observe from 
that curve that the attraction starts at about 1.7 fm instead of the 2 fm in the case of the total potential discussed 
above. This result is in some way in agreement with the phcnomcnological isoscalar spin-orbit potential for which it 



is well known that scalar- isoscalar and u> mesons are the mesons which play the most important role 1 1 1 . 

Similar results have already been found in Ref. fTlj] where a different way has been used to introduce a scalar- 
isoscalar degrees of freedom. In Ref. fl7|| , the scalar-isoscalar meson field is a dilaton, and it explicitly couples to the 
sixth-order term. However the contribution to the spin-orbit force coming from the sixth-order term was not complete 
in Ref. p7| . Indeed, these authors ignored the baryon exchange current contribution to the spin-orbit force in their 
calculations, and we have shown in Ref. that this contribution is significant, with respect to the direct term fll2] , |l7| ] , 
and should be taken into consideration. It is then legitimate to question whether the qualitative results of Ref. [17] 
remain valid if the baryon exchange current contribution is taken into account. In this work, by using a different way 
of including the scalar-isoscalar meson in the Skyrme Lagrangian [14,16], and by considering the entire contribution 
of the baryon current to the spin-orbit channel, we show that the main result of Ref. JTtJ (i.e., an attractive N-N 
spin-orbit force at large distance from Skyrmions with scalar mesons) is still valid. In this sense, my results can be 
viewed as a confirmation of those of Kalbermann and Eisenberg p"7| ] . 

Even though we have obtained an attractive force at large distances in the isoscalar spin-orbit channel, attraction 
is still missing at intermediate distances, and thus the problem remains unsolved within Skyrme-type models. This is 
not surprising since we have considered here only a finite-mass scalar field without the other vector mesons. However, 
the result obtained here, namely, the change from a zero contribution to the isospin-independent spin-orbit interaction 
to an attractive one when replacing the frozen e- meson field (nonlinear a model) with a realistic one (linear a model) , 
is very encouraging. Indeed, it suggests that the right way to obtain an attractive isospin-independent NN spin- 
orbit force within Skyrme-type models in the framework of the product ansatz is to replace the pion theory (|l|) 
with a realistic effective model including, in addition to the scalar-isoscalar e field, all low-lying vector mesons and 
taking into account the finiteness of their masses. In a sense, the model considered here can be seen as a minimal 
and modest improvement of the pion Skyrme theory. Similarly, it has been shown in Ref. ]l4| that the model (|^) 
gives rise to attraction in the NN central potential but the attraction occurs at distances larger than required by 



phenomenology. The problem has been finally solved when the other vectors mesons were included J_4 16 1. For this 
reason, we believe that in order to cure the problem of the isoscalar spin-orbit force, one has to consider effective 
models which incorporate the first mesonic resonances with finite masses. For instance, when considering a finite-mass 
w-meson model, the w-field couples directly to the nucleon via the baryon current (defined after Eq. ([!])) and generates 
the common three-pion exchange diagram, contrary to the case of the local approximation Cq (cf. Eq. (|l|)). Indeed, 
in the latter, an unexpected two-pion exchange piece coming from the baryon exchange current arises in addition 
to the three-pion piece and contributes with a positive sign to the isoscalar spin-orbit force yielding a repulsive 
interaction fl3|| . This problem with the baryon exchange current is obviously avoided when the sixth order term C§ 
is replaced with a more realistic w-meson model. 

Finally, in addition to considering finite-mass mesons Lagrangians, we would like to mention a further way that 
might lead to the desired attractive isoscalar spin-orbit force within the framework of Skyrme-type models. It is 
concerned with the approximation of the product ansatz configuration (1q). In the latter, the two nucleons are 
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supposed to keep their spherical shape without deformation even after overlapping. This is certainly not true in 
reality. Thus, while still using the product ansatz so that we keep benefiting from its simplicity and being able to 
perform analytical calculations, we may improve on the approximation (|^) by allowing the shape of the nucleon to 
deform when approaching the other nucleon (c/. Ref. j23j and references therein). Calculations of the spin-orbit force 
within this approach are under way p4j . 
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FIG.l. The chiral function F (full line) and the £-field (dashed line) in fm 1 solutions of Eq. (|j). The 
dotted line corresponds to the chiral field F in the nonlinear model (f = f n ). 




FIG. 2. The potentials Vlo- (dashed line), V4 (dotted line), Vq (dashed-dotted line), Vl ct + Vg (starry line) 
and Vl ct + V4 + Ve (full line) with respect to the relative NN distance r. 
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